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ABSTRACT 

A unified, three-scale system of equations accommodating nonlinear 
velocity driven helical dynamos, as well as time-dependent relaxation of 
magnetically dominated unihelical or bihelical systems is derived and solved 
herein. When opposite magnetic helicities of equal magnitude are injected 
on the intermediate and small scales, the large scale magnetic helicity grows 
kinematically (independent of the magnetic Reynolds number) to equal that 
on the intermediate scale. For both free and driven relaxation large scale fields 
are rapidly produced. Subsequently, a dissipation-limited dynamo, driven by 
growth of small scale kinetic helicity, further amplifies the large scale field. The 
results are important for astrophysical coronae fed with bihelical structures by 
dynamos in their host rotators. The large scale for the rotator corresponds to 
the intermediate scale for the corona. That bihelical magnetic relaxation can 
produce global scale fields may help to explain the formation of astrophysical 
coronal holes and magnetohydrodynamic outflows. 

PACS codes: 52.30.Cv, 95.30. Qd, 52.65Kj, 96.60.Pb, 98.35.Eg, 98.62Mx 



I. Introduction 
A. Background 

Large scale helical dynamo theory provides a framework for modeling the in situ origin 
of large-scale magnetic field growth in planets, stars, and galaxies^"^ and has also been 
invoked and to explain the sustenance of fields in fusion devices. Much of the associated 
dynamo research can be divided into (1) that which parameterizes the nonlinearities in 
an effort to practically model the field evolution of specific objects or systems, and (2) 
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that aimed at identifying the physical principles of the nonlinear dynamo by dynamically 
including the mutual backreaction between the field and velocities for simplified systems. 

Progress in understanding the nonlinear helical field growth and saturation has resulted 
from a combination of numerical and analytic work that dynamically incorporates magnetic 
hehcity evolution. ^^"^^ Analytic models^^"^^ employing a discrete number of scales and 
including the backreaction of the field growth on the turbulent velocity have shown that 
for kinetic helicity driven dynamos, the build up of small scale magnetic helicity ultimately 
quenches the large scale field growth, in agreement with numerical simulations. A 
two-scale nonlinear theory of large scale field growth from a dynamo driven by small 
scale magnetic helicity fluctuations rather than kinetic helicity fluctuations has also been 
developed. In this case the kinetic helicity becomes the quenching agent rather than the 
growth driver. Like the velocity-driven helical dynamo, such magnetically driven helical 
dynamos thrive on a finite turbulent electromotive force in Ohm's law. 

A magnetically driven helical dynamo can also be described as "dynamical magnetic 
relaxation." Magnetic relaxation describes the process by which magnetic structures in 
magnetic pressure dominated environments evolve to their equilibrium states. The fully 
relaxed state is the Taylor state, determined by minimizing the magnetic energy subject 
to the constraint that magnetic helicity is conserved. The result is a force-free helical 
configuration with the field relaxing to the largest scale available, subject to boundary 
conditions. However, Taylor's theory by itself is not a dynamical theory since it does not 
provide a time-dependent description of how the large scale magnetic helicity evolves. A 
time-dependent theory that also includes velocity and fluctuations is required. 

In fusion devices, where the plasma is magnetically dominated, the word "dynamo" 
is typically used to describe a magnetically dominated dynamo in which magnetic flux 
is converted from toroidal to poloidal (or vice versa), increase the scale of the flcld, and 
sustain a strong magnetic flux against microphysical dissipation. ^^^^ Detailed studies of 
these dynamos focus on the instabilities that drive the fluctuations (tearing modes) as well 
as the subsequent relaxation, conversion and sustenance of the magnetic flux. The full 
boundary value problem can be solved for speciflc conflgurations. 

Similar magnetically driven dynamos are likely an important part of the generation 
of large scale astrophysical flelds that mediate coronal holes and magnetohydrodynamic 
(MHD) jets in astrophysics. The analogy can be drawn between spheromak or reversed-fleld 
pinch formation^^ and coronal fleld relaxation: The laboratory conflgurations form because 

the injection of magnetic helicity into a system takes the system away from the relaxed 
state, but also drives tearing mode fluctuations that yield a flnitc turbulent electromotive 
force. This in turn drives the system to to a new relaxed state. In astrophysical coronae 
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the situation is similar. The helical field produced by a velocity driven dynamo inside of a 
rotator is buoyantly fed from below into the corona. This takes the magnetically dominated 
corona away from the relaxed state initially, but fluctuations arise in which the system is 
then driven back to a relaxed state. The corona therefore acts as the "laboratory" plasma 
volume, injected with helicity from below. In the astrophysical case, the injection occurs 
not at one site, but throughout the coronal-surface boundary. While each of the individual 
injection sites (flux tubes) would be analogous to a single laboratory conflguration if the 
twists are injected along the tube, each tube can also open up to even larger scales by 
mutual interaction. In this way, magnetic energy can appear on scales even larger than 
those on which the fleld was originally present. This is a difference between the laboratory 
and astrophysical cases. 



B. Astrophysical motivation for bihelical relaxation 

Magnetically driven dynamos are important for a wide range of astrophysical coronae 
because all coronae of turbulent astrophysical rotators (the sun, stars, and accretion disks) 
are likely magnetically dominated. For the solar corona,^^ observations support the basic 
paradigm that particle acceleration and heating result from extraction of coronal magnetic 
energy via reconnection or wave dissipation. The coronal magnetic fields come from 
magnetic loops produced within the sun that are injected into the corona. Because the field 
polarity reverses every 11 years, these loops could not have been merely the residue of the 
accreted interstellar field. The sun also has much larger scale ("open") fields in coronal 
holes along which the solar wind propagates. These fields also reverse during the solar 
cycle, again highlighting that they are not simply the result of flux freezing from the sun's 
formation. That the geometry and topology of smaller bipolar structures evolve to produce 
global structures in the solar cycle can be modeled empirically.^^ 

Accretion disk coronae likely involve similar processes. Non-thermal emission and 
flares from accreting stars or black holes are also naturally explained by the emergence of 
magnetic flelds from the turbulent plasma below.^^"^^ By analogy to the sun, disk flelds 
supplied to its corona can open up. Because accretion disks arc "fast rotators" (as are young 
X-ray luminous stars'^'') the open field lines are not only the track along which outfiows 
propagate, but can act as a drive belt to extract the rotational energy in the disk (or young 
star) supplied by accretion. Conversion of the accretion energy fiux into Poynting fiux can 
drive the observed outfiows. ^''"^^ Though particular MHD models differ as to how far above 
the disk the outflows are carried as Poynting flux,^°' they all involve large scale flelds and 
collectively offer the most viable jet/outflow paradigm for a range of accretion engines. 
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MHD jet/outflow studies typically impose the global field anchored in the turbulent 
rotator as an initial condition, and then allow the system to evolve or consider the steady 
state. Because field reversal cycles cannot be directly measured in accretion disks, 
it is not as straightforward as it is for the sun to directly determine whether the large 
scale fields are accreted or produced in situ. But difficulties in accreting ffux in a thin 
turbulent disk^^'^^ and the fact that the solar corona produces large scale global fields in 
situ provide reasons to believe that in situ production in disk coronae is ongoing. The 
supply of magnetic energy to a corona is directly supported by 3-D MHD simulations,^^ but 
the self-consistent generation of large scale flux still needs further study. 

In the standard use of helical dynamos to generate large scale astrophysical fields, an 
initially weak field is assumed and the growth of the large scale field is studied for given a 
choice of velocity driving and global boundary conditions.^'^^'^^'^^"^^ While this approach 
applies inside the rotator, it hides an important aspect of the large scale field generation 
in the corona. The large scale coronal fields might be more appropriately described as 
the second phase of a two-phase process: (1) In the first phase, a velocity driven helical 
dynamo generates large scale fields inside the rotator. The field produced then emerges into 
the corona. (In order for buoyancy to beat shredding for a disk, a helical dynamo within 
the rotator seems to be required because non-helical turbulent amplification produces 
structures that are dissipated before they can rise significantly'^'^.) (2) Because the "large" 
scale field produced by the velocity driven dynamo in the disk is of "intermediate" scale 
with respect to the corona, it must incur further relaxation to open up to the global scales 
required by holes and outfiows-just as we see in the sun. However, this relaxation occurs 
in an already magnetically dominated environment, so phase (2), just discussed, is indeed a 
magnetically dominated dynamo. Because magnetic helicity conservation implies that any 
nonlinear dynamo in the rotator below produces helical fields of opposite signs on different 
scales, i5~23 j^ggg pjg^ ^ Ylef. 23 for an illustration) the question becomes: Is a net global 
large scale field produced when injected with bihelical fields that supply a zero net helicity? 

C. Present aim: closed time-dependent theory, not detailed boundary value 

study 

Unlike the laboratory dynamo, the astrophysical coronal dynamo involves multiple 
helicity injection sites (i.e. each buoyant twisted fiux tube) rather than a single injection 
(helicity injection along the mean field.). Nevertheless, in a real corona, like in a laboratory 
configuration, the global magnetic relaxation will result from current instabilities instigated 
by the magnetic hehcity injection. Indeed, a full treatment of coronal bihelical magnetic 
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relaxation will ultimately require boundary value studies. However, much insight into 
the basic nonlinear effects of magnetic helicity transfer and velocity backreaction can be 
gained using a minimalist mean field approach without a detailed modeling of the tearing 
modes that drive the fluctuations. Rather than a detailed boundary value problem, I 
consider a simpler problem in which the scale magnetic helicity is driven into the system 
quasi-isotropically at the intermediate and small scales and the subsequent growth at 
at the large scale calculated. This is directly analogous to the approach used to study 
nonlinear principles of velocity driven helical dynamos, and has led to significant progress in 
reconciling results from analytic treatments and numerical simulations.^^' Keeping 
the set-up simple (a closed system driven quasi-isotropically with kinetic or magnetic 
helicity fluctuations), facilitates calculation of how much relaxation occurs on dynamical 
vs. dissipativc time scales, the strength of the backreaction on the relaxation from velocity 
fluctuation growth, and the properties of the saturated state. 

Toward this end, I derive and solve a set of equations that dynamically models the 
time evolution of a system injected with bihelicity (opposite signs of magnetic helicities on 
different scales) to determine the extent of the magnetic helicity generation on even larger 
scales. The system is magnetically dominated in the sense that the energy injected into the 
system is initially entirely magnetic, and the total magnetic energy (summed over all scales) 
remains dominant throughout the calculations even when the velocity fluctuations grow. 
The time dependent evolution requires a three-scale nonlinear theory. The nonlinearity 
arises because the magnetic helicity equations are coupled to dynamically to equations for 
the velocity fluctuations. The latter are generated by the magnetically dominated system 
on its journey to the relaxed state. 

In Sec. II, I derive the equations and collect them in dimensionless form. In Sec. Ill, 
I solve the equations and discuss the solutions for the following cases: (1) driven bihehcal 
magnetic relaxation, where the system is continuously driven with oppositely signed 
magnetic helicities on different scales; (2) kinetic helicity driven dynamo, where the system 
is continuously driven with kinetic helicity only on the intermediate scale; (3) unihelical 
magnetic relaxation, where the system is continuously driven with 1 sign of magnetic 
helicity only on the intermediate scale; (4) free bihelical magnetic relaxation, where the 
system is initially injected with oppositely signed magnetic helicities on the intermediate 
and small scales but with no subsequent driving. In Sec. IV, I summarize the results and 
the implications. 
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II. Equations for bihelical magnetic relcLxation 

In this section, the generahzed set of equations that can describe both kinetic hehcity 
and magnetic hehcity driven dynamos is derived. 



A. Magnetic helicity equations 

I assume that there are 3 scales on which magnetic hehcity can reside, defined such 
that their wave numbers satisfy > k2 > ki > 0. I take B, A, and E to be the components 
of the magnetic field B, vector potential A, and electric field E that vary on the scale ki^. 
The respective quantities on scales ^ and k^^ will be labeled with subscripts 2 and 3 
respectively. 

In a closed or global volume, the total magnetic helicity, H — ( A • B) , satisfies^'^^'^^ 

dtH ^-2{E-B), (1) 

where E = —dtA — V0, 4> is the scalar potential, and () indicates a global volume 
average. Such a conservation formula also applies separately to each scale. That is, using 
E = —dfA — V0, 62 = —dtai2 — V02 and = — S^aa — V03 and dotting with B, b2 and ba 
respectively we obtain 

dtH, = -2(E-B) (2) 
dtH2 = -2(e2 • b2) (3) 

and 

dtHs = -2{es-hs), (4) 

where Hi = {A ■ B), H2 = (a2 • b2), and H3 = (as • hs). The expressions for (E • B), 
(e2 • b2), and (es • ba) are 



(E-B) = -(£:-B)-(e-B) + A(J-B) (5) 

(e2 • b2) = • B) - (£2 • b2) + A(j2 • b2) (6) 

and 

(es • bs) = (e2 • b2) + (e • B) + A(js • bs) , (7) 

where A is the magnetic diffusivity, J = VxB, j2 = Vxb2 and js = Vxbs are the 
normalized currents on the respective scales defined in convenient units, £ = V2 x b2, 
e = V3 X bs, €2 = (vs X b3)2, and ()2 indicates averages which vary on scales > /c^^ and V2 
and Vs are the associated velocities. 
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Combining equations (2-4) with (5-7) gives 

dtHi = 2E\\\Hl^^\ky^ + 2e\\\Hl^^\kY^ - 2XklHi (8) 
dtH2 = -2?|| l/f^'l A;^' + 2€2iiH^^^\kl^^ - 2XklH2 (9) 

and 

dtH^ = -2e\\\Hl'^\k\'^ - 2e2\\Hl'^\k\'^ ~ 2\klH,,, (10) 

where S\\= £ ■ -jSj-, e|| = e ■ , and e2\\ = 62 ■ Here the three-scale approach has been 
used to write (J ■ B) = kfHi, {^2 ■ = ^2-^2, and (js ■ h^) = k^H^, and maximal helicity 
has been assumed so that (B ) = kjHi, (b|) = A;|i72, (bg) = k^Hs. 



B. Time-dependent electromotive force equations 

A dynamical nonlinear helical dynamo theory should include^^'^'' the time evolution of 
£ and by generalization to three scales, also the time evolution for e, and €2- Of primary 
interest are the specific components S\\, ey and e2\\. We therefore have 

B B B 

dtS\\ = {dtV2 ^ t>2) • PI + (v2 X ^tba) • ^ + £ • ^t^, (11) 

B B B 

dte\\ = {dtVs X bg) • ^ + (vg X dths) • ^ + e ■ dt-=- (12) 



and 



5*6211 = {dtVs X b3)2 ■ + (V3 X dth3)2 ■ ^ + €2 • (13) 



To proceed, we need equations for 9tb2, dfhs, dtV2: and dtVs, which come from the magnetic 
induction equation 

a^B = V X (v X B) + AV^B, (14) 
and the incompressible Navier-Stokes equation 

a^v = v X u; - V(P + vV2) + J X B + i/VV, (15) 

where P is the thermal pressure, and a; = Vxv. Because there are multiple fluctuation 
scales {k2^ and ^3"^), we must perform multiple averaging. To obtain the equations for 
quantities with subscript 2 and 3 note that 

dth2 + dth3 = {dtB-dtB), (16) 
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so that 

dth2 = {dtB - dtB)2, (17) 

and 

dths = {dtB - dtB) - {dtB - dtB)2. (18) 

The analogous relations apply for dtV2 and dtVs, with the simphfication that dtV — for 
the present calculation. Thus 

9tV2 + dtV3 = dtv, (19) 

so that 

dtV2 = {dtv)2, (20) 

and 

dtVs = dtv - {dtv)2. (21) 
After some algebra, applying this procedure to (14) and (15) gives 

dth2 = VX(V2 X B) + VX(V2 X b2)2 + VX(V3 X b3)2 + AV^b2 + f62, (22) 



dths = VX(V2 X ba) + VX(V3 X b2) + VX(V3 X B) + VX(V3 X b3)3 + XV^hs + fb3, (23) 

{^t^r2)^ = n,,(B • Vb2 + b2 • VB),- + u{V\2)i + f.2, (24) 

and 

{dtV3)i = Uij{B ■ Vb3 + b3 • VB + b2 • Vb3 + b3 • Vb2),- + iy{V\3)i + f,3, (25) 

where u is the viscosity, ffc2 and f62, are forcing functions assumed not to correlate with 
Vi,V2, and fy2 and £^3 are forcing functions assumed not to correlate with b2, b3. The 
projection operator Iljj = (5ij — ^'Y^ ) comes from taking the divergence of the velocity 



equation to solve for the total pressure of incompressible flow. 

In deriving (17-25) and in subsequent manipulations, the following important 
simplifying assumptions are made: (1) Mixed correlations between quantities of widely 
separated scales arc assumed small. Thus I ignore terms of the form e.g. (j2 ■ h^). (Note 
however that e.g. 32 ■ bs 7^ 0; only when put in brackets is the assumption applied.) (2) 
The magnetic field or velocities on a given scale are assumed to be fully hehcal, so that 
specifically V3 x u?3 ~ V2 x a;2 — V x (VxV) ~ x bs ~ j2 x b2 ~ J x B ~ 0. This 
is not as restrictive as it might appear because the mixed quantities e.g. j2 x are not 
assumed to vanish and thus non-linear couplings and transfer between scales are included. 



-9- 



(3) Cross helicities are ignored (e.g. pseudo-scalar mixed correlations between v and b such 
as (b • v)). (4) It is assumed that terms involving the correlation of products of separate 
quantities averaged over the intermediate scale can be factored into two global correlations. 
That is, ((ja • bsja^i) ^ {h ■ hs){bl) and ((vg • u;3)2&i) (vg • uJs){bl). 

Armed with these assumptions, Eqs. (22-25) must now be used in Eqs. (11) and (13). 
Eqs. (22) and (23) are required for the third terms in (11) and (13). These terms can be 
computed in configuration space by direct substitution of the latter two equations into the 
former. Eqs. (24) and (25) arc required for the second terms in (11) and (13). These terms 
are computed by Fourier transforms, and triple correlations are also retained by use of 
the "minimal r" closure '^^ (which obviates the first order smoothing approximation) by 
approximating the triple correlations as the product of the respective electromotive force 
divided by a damping time scale. Applying this closure here, the triple correlation arising 
from the fifth term of (23), when input into the third term of (13), contributes a term of 
the form — enlvs/csl, and the third and fourth terms of (22) when plugged into the third 
term of (11) contribute terms of the form --S\\\v2k2\, and — e|||v3/c3| respectively. It should 
also be noted that the last terms of (11) and (13) respectively will also be included in these 
damping terms. The coefficients of the damping terms can in principle be taken to be 
different from unity^° but this is not considered further here. 

Using the above approximations and plowing through the algebra that arises from 
plugging Eqs. (22-25) into (11-13) gives 

dt£\\ = (klH^ - H^) f -^^k - ^11^2(^2/^2 + A + i/) - £2,1^2^2, (26) 



dfe, = (klH, - H^,) ^ilJ^ - - e\\kl{v,/k, + A + .), (27) 



and 

|rTl/2|,l/2 ttV t.S/'i TT 

9,6,1 = {klH, - H^) - - e2,kl{v,/k, + A + .), (28) 

where H2 = (v2 • UJ2) and = (V3 • a;3) are the kinetic hehcities associated with the 
respective scales. 



C. Kinetic helicity equations 

To close the set of equations, we need dynamical equations for H2 and H^. Using (20) 
and (21) in (15) to obtain equations for dtV2, dtU}2, dfVs and dtVs, ignoring total divergences 
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inside of averaging brackets and using the assumptions discussed below (25) we obtain 

dtH^ = 2(J ■ (b2 X a;2) + 2(B ■ (wa x h)) ^ 2(A;i - kik2){S ■ B) - 2kluH^ (29) 

and 

dtH^ = 2(3 ■ (ba X a;3) + 2(B ■ {cv^ x j3)^+ 2(j2 ■ (bg x a;3)) + 2(b2 ■ (c^s x js)) . . 

~ 2{kl - kik^){e ■ B) + 2{kl - k2ks){e2n ■ b2) - 2kli^H^ , ^ ^ 

where the latter similarities in (29) and (30) follow from maximally helical assumptions 
j2 oc b2, oc V2, js oc h^, and u:^, oc V3, and the assumptions that = ^j]^ and 

= r^rfeH- These latter relations will prove to be self- consistent for the solutions 
discussed later. 



D. Aggregate set of equations in dimensionless form to be solved 

Before collecting the complete set of equations to be solved, it is useful to add injection 
and loss terms to the three magnetic helicity equations derived above. Since we will study 
the evolution of magnetic helicity for a system subject to injection of bihelicity — that is, 
helicity injected of opposite sign on the scales k2^ and k^^ — we add terms to the H2 and 
H3 equations that account for injection (or loss) of magnetic hehcity on these respective 
scales. In addition, if the system is a corona of finite volume, large scale fields therein might 
be able to escape as Poynting flux. While the injection terms just described can be used 
to combine injection and loss for H2 and H3, for Hi I add a distinct Poynting loss term 
oc V^i^i/ci ~ H^^'^kl^^, where Ya = B is the Alfven speed associated with B (recall that 
magnetic fields are written in Alfven units herein). (When losses are included, the magnetic 
helicities should be re-interpreted as gauge invariant relative helicities^^ with respect to a 
background potential field, but the basic formalism herein does not change significantly.) 

For completeness, I also include injection (or loss) terms in the H2 and equations. 
Use of the injection term for H2 will help illustrate the relation of the present study to the 
2-scale solution of the kinetic helicity driven dynamo.^'' 

To write the needed equations in dimensionless form, define hi = Hi/\H2fi\, 
h2 = H2/\H2,o\, = Hs/\H2,o\, hX_= /\klH2fil = /\klH2fll Rm = Hl^^lXkl'^ 
Rv = hIJo / ukl^"^ , T = tkl/2H2fl, Q = £^||//c2i?2,o, q = e||/^_2^^2,o, and q2 = e2||//i;2^^2,o, where 
H2fi = -^2(0) if magnetic helicity is initially injected and i72,o = -f^2'(0)/^2 if kinetic helicity 
is initially injected. Using these dimensionless parameters and adding the afore mentioned 
injection and loss terms, Eqs. (8-10) become 

= 2Q|hl'^| {^f + nhl"\ (fe)"' - £ - c:|fcf (31) 



drh, ^s,- 2Q\hr\ {^y' + 2q,\hy'\ - (32) 

and 

9M = .3 - 2u\hr\ [^f - 2,,\hr\ - t; (U)' , (33) 
where ci is a constant parameterizing the Poynting loss term and Si and S2 are the hehcity 
injection terms described in the previous paragraph. 

Completing the set of equations using the dimensionless variables defined above, we 
have for Eqs. (26-28) and (29-30) 

(34) 




5r?2 = -(/i3 - hl)\hy'^ 




- -\h 




dX^sy^2Q\hr\{^f\l-^^y% (37) 

and 

B.hX = 4 + 2 (I - if) (Ij)"^ + 2 (I - I) - % [^f , (38) 

where and s\ are the analogous source terms to S2 and S3. Eqs. (31-38) can be solved 
for /ii, /i2, ^3, Q-i Q, 0.2-1 , and as a function of time for a variety of different initial 
conditions. A few cases are studied below. 



III. Discussion of solutions 



Eqs. (31-38) provide a unifying framework for both velocity driven and magnetically 
driven nonlinear dynamos. They characterize the magnetic and kinetic helicity dynamics 
for a range of different physical scenarios, depending on the initial conditions. In all cases, 
below ki — 1, k2 — 5 and ks — 20. 
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A. Driven dynamical bihelical magnetic relcixation 

Consider the case for which dth^ = dth2 = s\ = s\ = and /i2 = — /is = 1 for all times. 
This implies continuous injection and removal of oppositely signed magnetic helicities on the 
intermediate and small scales so as to maintain h2 = —h^. Fig. 1 shows the solutions of Eqs. 
(31-38) for initial {t = 0) values hi{0) = 0.01, h^{0) = h^{0) = Q{0) = q{0) = ^2(0) = 0. 
Fig la shows hi through the saturated regime and Fig lb focuses on its early time solution. 
Fig Ic shows /i2 and , Fig Id shows hs and h^, Fig.le shows Q,q, and 52, and Fig If is a 
diagnostic solution. These are discussed below. 



1. Late time regime 

To analytically estimate the late time steady state solution for hi of Fig la, note that 
from equation (37) in the steady state, Q is inversely proportional to the large quantities 
Rv and hY^. Assuming that \Q\ « \q\ and \Q\ « \q2\ (these will turn out to be vahd 
assumptions), we can then use Eqs. (31), (36), and (38) to find the saturation value of hi. 

In the steady-state, Eq. (38) gives 

assuming that the last term in (38) is small. The fact that q2 asymptotes to negative values 
imphes that hs ~ h^ in (36) due to the growth of h^ ■ The steady hmit of Eq. (36) can 
then be rearranged to give the approximate relation 

* - - jf. («) 

where we have neglected the Rm and Ry terms in that equation as they are small. 
Combining (39) and (40) gives 

7.2 \ 1/2 /u \ 1/2 /-, ,^ \ 

(41) 

This equation can then be used in the steady-state limit of Eq. (31). Using the earlier 
assumption that \Q\ << \q\, we then have 

^, Rm kj ( l-k2/h \ 

3 klkl \l-ki/k,)' ^^'^ 
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which provides an excellent approximation to the actual solution at late times. It also 
shows the insensitivity to Ry- 

The saturated hi oc Rm in (42), is also exhibited in Fig. la. This can be understood 
as follows: hi grows first with the sign of /i2 by direct analogy to unihelical magnetic 
relaxation/^ (see also sec. 3.2.2): The magnetic helicity seeks the largest scale available of 
a given sign. However in the bihclical case, the additional injection of negative drives the 
growth of negative /i^. In fact, \h'^\ grows to very slightly exceed {h^l as seen in Fig. Id. 
Then, q becomes slightly positive. This drives additional hi growth by an induced kinetic 
helicity driven dynamo from h^ , by maintaining the q term on the right of (31) at a positive 
value (long dashed curve in fig le). The larger the value of Rm, the less competitive the 
last term on the right of (31) is with the q growth term. This explains the Rm dependence 
of the saturation value of hi. 

The above interpretation can be tested by artificially reducing the coupling of h^ to 
the other equations. This was done to produce the solutions in Fig. If by multiplying the 
entire right side of (38) by 10~^. Notice the significant reduction in the value of hi at even 
relatively early times compared to Fig la. In this case, the presence of h^ acts to damp the 
growth of hi because q remains negative at all times. The two curves in Fig le also confirm 
that it is indeed the 52 term in (34) that leads to the rapid growth of hi in Fig la,b since 
lowering the coefficient on that term lowers the saturated hi. 

2. Early time, kinematic regime 

In Fig. lb, the early time solution for hi is shown. The solution is independent of Rm 
up to /ii ~ 1. At early times, growth of hi is dominated by the first term on the right side 
of (31). The dominant term for the growth of h^ from (37) is very similar and thus the two 
initially grow at nearly the same rate. However, h^ eventually builds up to a value that 
stops the growth of Q in (34), after which the latter decays and the growth of hi and h\ 
slows. The decay of Q occurs when the first two terms on the right of (34) approximately 
balance. Since /12 = 1, the critical /i^ is ~ /i2/ (1 + ki/k2) ~ 0.83. This provides an excellent 
match to the value of hX at the end of the kinematic regime shown in Fig le. Since the 
growth rate of h^ is slightly less than that of hi due the presence of the (1 — A;i/A;2) factor in 
the first term on the right of (37) compared to the first term on the right of (31), hi grows 
a bit larger than /i2 ~ 1 during the kinematic regime. After this time the growth depends 
on Rm. Just after the kinematic regime, the hi curve for lower Ry is the higher of the two 
curves in Fig. lb. This is because a lower Ry reduces the kinetic helicity backreaction in 
Q since the velocity decays faster. However, the two curves eventually cross to produce the 
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late time solution of Fig. la for the reasons discussed in the previous subsection. 

The value of hi reached during the kinematic regime during bihelical relaxation is 
approximately the same as that reached in the kinematic regimes of the simpler two-scale 
unihelical magnetic relaxation modeP^ and the kinetic helicity driven helical dynamo^° for 
the similar reason that in these cases, growth of the largest scale hehcity is driven by the 
difference between magnetic and kinetic helicities from the small scale. Driving a two scale 
system with kinetic (magnetic) or magnetic helicity grows the large scale hi rapidly at 
first, but then the magnetic (kinetic) helicity on the small scale builds up and quenches the 
growth to dissipation limited rates. As discussed in the previous subsection, the evolution 
of hi differs for bihelical relaxation long after the kinematic regime. Sec. III.B shows how 
all of these cases are consistent and just depend on the choice of initial conditions. 

Finally, a few more comments on Fig. Ic and d: Although curves for both Rm = 1000 
and Rm = 6000 are shown in both panels, they overlap at early times when there is no 
Rm dependence. Also notice the oscillations in h^ , which damp very quickly, after which 
h^ remains steady. Remember that the unit of time is defined on the A;^^ scale. The short 
evolution time of h^ results from the factors of {ks/k2y — 16 in Eq. (38). 

3. Comparing k2 — vs k^ > k2 

If k2 = ks, then for the case of Fig. 1 where /i2 = — /is = 1, the net injection of helicity 
at the scale k2 — k^ would be zero and nothing should happen. This has been tested and 
indeed there is no growth for this case. This highhghts that although growth of hi can 
occur when the net injection of h2 + hs — 0, the two signs of helicity must be injected on 
sufficiently different scales (A;|/A;| >> 1) for non-trivial solutions. 

4- Inclusion of boundary terms 

In Fig.l, ci = 0. Fig. 2 shows the effect of varying ci > 0, to approximate different 
finite loss rates of hi . The result is that the value of hi saturates at much lower values than 
in Fig. la, as expected. 
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5. Case of fixed dth2 — —dfh^ — constant 

In the previous subsections, S2 = S3 = and h2 = 1 and h-^ = — 1 were held constant. 
The case of constant injection rate, namely Si = 1 and ^2 = — 1, is a natural alternative to 
consider. Solutions for this case also reach a steady state at late times. Higher values of hi 
are obtained in this case than for the fixed h2 — —hs case. I do not discuss this case further 
here as the implications are similar. 

B. Kinetic helicity driven dynamo and driven unihelical relcixation 

In previous two-scale and three-scale approaches for the kinetic helicity driven 
dynamo, ^''"^^ dynamical equations for the kinetic helicity were unnecessary because the 
kinetic helicity was kept steady. In contrast, the unihelical relaxation study did include 
kinetic helicity dynamics but it was a two-scale approach with only 1 sign of injected 
magnetic helicity. That the equations (31-38) are more general means that they should 
reduce to these other cases in the appropriate limits. 

1. Kinetic helicity driven dynamo 

If the system is driven steadily with h^ starting with only a seed hi and /i2 = /13 — 0, 
then the solution reproduces the non-linear multi-scale helical dynamo^^, which in turn is 
consistent with the results of the two-scale version. ^° Fig. 3 shows the resulting solution 
for hi, which represents a simple dynamo. The rate of growth up to /ii = 1 is fast 
and independent of Rm- The ultimate steady state ratio /ii//i2 = {k2/kiY = 25 to which 
the solution converges at late times [t > -Rm^2/^i), is also independent of Rm- This can be 
seen from solving equations (31) and (32) with S2 = ci = g = ^2 = in the steady state). 

The rate of approach to the steady state is slower for smaller Rm- This is because the 

growth of hi for the /i^ driven dynamo occurs by segregating opposite signs of magnetic 
helicity to small and large scales. Growing hi larger than /i2 means that h2 has to dissipate 
faster than hi- The larger the Rm, the slower the dissipation of /i2, leading to a smaller Q, 
and a longer time for hi to saturate. 

This late-time dependence on Rm differs from that of hi growth in bihelical magnetic 
relaxation as discussed in Sec. III.A.l. There, a larger Rm gave a larger saturation value 
and a faster growth rate at late times. The difference arises because growing hi at late times 
in that case is not determined by losing small scale magnetic helicity but by maximizing 
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the effect of the secondary induced driven dynamo. There is no contradiction in the role 
of Rm between these two cases because of the difference in the initial conditions. 

The saturated value of hi depends slightly on the value of turbulent diffusion (the 
third term of (34)) and its quenching. Different scaling factors for turbulent diffusion were 
used previously^"' and these differences explain the small difference in the saturation 
values therein when compared to Fig 3. There is no dynamical quenching of turbulent 
diffusion in either the present or previous approaches, so here I simply eliminated the extra 
multiplicative phenomenological factors. 

In short, the equations of Sec. 2 fully account for results in the previous studies of the 
nonlinear dynamo. 

2. Unihelical magnetic relaxation 

If the system is driven with /i2 but without hz, then the growth of hi responds just 
as in unihelical dynamical relaxation. The injected helicity /i2 drives the growth of hi by 
dynamical relaxation. (This is called unihelical relaxation because the helicity is injected 
with 1 sign at h2 and the same sign then grows at hi ) Fig. 4 shows the growth of hi. 
The kinematic regime proceeds until hi ~ 1, after which the growth of hi depends on 
dissipation. However, unlike the bihelical magnetic relaxation of section 3.1, here the 
late time saturation value of hi depends on the ratio Ry/ Rm- Here /i^ plays the role of 
back reactor, and a larger Ry means that hX is more effectively drained, allowing further 
hi growth until saturation. Again note that in the contrasting bihehcal case of Fig. 1, 
the saturation value of hi is insensitive to Ry (for Ry >> 1). There is no contradiction 
because in that case, both positive and negative signs of magnetic helicity are injected and 
the h^ = —1 injection produces an additional h^ driven dynamo through q in (31) that 
dominates any effect of the suppression of Q at late times. 

The important point is that the previously studied case of two-scale unihelical 
relaxation emerges from the more general equations of Sec. 2 by the appropriate choice of 
initial conditions. 



C. Free dynamical bihelical magnetic relaxation 

In the previous sections, either kinetic or magnetic helicity was assumed to be 
continuously driven into the system. Here I consider the case where with initial 
^2(0) = — /i3(0) = 1 and hi{0) — 0.01, but with no subsequent driving or boundary terms. 
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i.e. Ci — S2 — — S2 — — 0. The system is then allowed to freely relax. Results of this 
case are shown in Fig. 5a-f. 

Prom Fig. 5a and 5b, it is evident that significant growth of hi occurs. The kinematic 
phase (independent of Rm and Ry) ends when /i2 ~ ~ ^1 — 1/2 by analogy to free 
unihelical relaxation/^ due to transfer of /i2 to hi. Fig. 5b shows that hi ~ 1/2 rather than 

/ii ~ 1 is reached in the kinematic phase because h2 also drives the growth of /i^. When 
~ 1/2, transfer of /i2 to hi is shut down as Q depletes from the decrease in the first 
term on the right of (34). After the kinematic phase, the h}^ induced secondary dynamo 
described in Sec. III. A takes over to amplify hi a bit beyond the kinematic value. Here this 
has only a small effect because without forcing the system eventually decays. 

For larger 7?^! the slower the eventual decay of hi. This is as expected from (31), but 
with the subtle implication that the Q and q terms combine to give a net contribution that 
does not affect the sensitivity to Rm from the dissipative term. For fixed Rm but smaller 
Rv, the decay rate of hi is the same, but the maximum hi reached is smaller. This is 
evident from comparing the middle and top curves of Fig. 5a and is correlated with the 
faster depletion of h^ and h^ in Fig. 5c and 5e for the lower Ry case. Remember that h^ 
drives a secondary dynamo that grows hi beyond the kinematic regime. As Ry is decreased, 
h^ decays faster making this secondary dynamo is less effective. Because /?.3, and thus 
are not continuously supplied, the secondary dynamo decays and so does hi. For arbitrarily 
larger Ry than shown in Fig. 5a, the maximum reached by hi before it starts to decay is 
hi - (1-A;2/A;3) =0.75< /i2(0). 

The curves for h-i, and h^ indistinguishably overlap in Fig. 5c but a blow up in Fig 
5d shows that the /i^ curves lie slightly to the right of the h^ curves. This slight offset 
is enough to drive the dissipation limited secondary dynamo just described. Fig. 5c also 
shows that h^ and /13 deplete faster for smaller Rm. This is expected from Eqs. (33) and 
(38) but with the implication that the q and terms combine to give a net contribution 
that does not contradict the Rm and Ry dependence from the dissipative terms. 

In contrast to the slower decay of hi, h^ and h\ with increased Rm, Fig. 5b. shows that 
/i2 and /i^ actually decay faster with increasing Rm- At first this may seem counterintuitive 
but it results because q2 is negative and decays more slowly with increasing Rm- Thus the 
contribution from the q2 term in (32) acts as a decay term. That this term is larger for 
larger Rm is evident from Fig. 5f where the bottom two curves show q2 for Rm — 2000 and 
Rm — 10000 respectively. The bottom two curves are close to the asymptotic values for 
large Rm- Little change occurs for higher value Rm > 10000 (not shown). 

Despite the subtleties, the key point to take away from the free bihelical relaxation 
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solution is that hi grows fast and decays more slowly than /i2 or h^. Like the driven case 
of Sec. III. A. when the system is injected with magnetic helicity of equal magnitude but 
opposite signs on intermediate and small scales /c^^ and k^^, significant magnetic helicity 
grows on the largest k^^ scale kinematically. For injected values of h2 = —h^ = 1, hi 
grows to hi ~ /i2(0)/2 ~ 1/2 kinematically and independent of hi{0). The growth time is 
approximately the Alfven crossing time of the fcf ^ scale using the injected field strength 
from /i2. Unlike the driven case where the maximum hi reached at late times is proportional 
to Rm, in the free case the maximum post- kinematic hi reached for Rm — > oo is /ii ^1. 
Subsequently, hi resistively decays. 

The free bihelical relaxation treated in this section is different from that of Ref. (52). 
There, the system was initially injected with helicity of one sign on the largest scale and 
the opposite sign on a smaller scale. In contrast, here the helicity is injected on two scales 
both smaller than the largest scale and hehcities on all three scales are allowed to evolve. 
Once hi grows to its maximum and h2 depletes, the subsequent resistive evolution of hi 
and hs acts qualitatively like the relaxation of the random bihelical case of Ref. (52) in that 
hi resistively decays more slowly than h^. 

IV. Summary and Implications 

A generalized set of equations that encompasses non-linear velocity driven dynamos 
as well as magnetically driven dynamos (equivalently, dynamical magnetic relaxation) is 
derived and solved herein. Previous work focused on these processes separately, but here 
the different processes are unified into a single set of equations that includes the dynamical 
evolution of kinetic helicity, magnetic helicity, and the turbulent electromotive forces. 
Particularly important is that the present approach includes the dynamics of 3 scales: 
two fluctuation scales in addition to the global scale, rather than just a single fluctuation 
scale and global scale. This is fruitful because it allows studying bihelical relaxation, the 
evolution of the helicities on all three scales when helicities of opposite signs are injected on 
the intermediate and small scale respectively. 

Example solutions of the bihelical magnetic relaxation equations reveal that even when 
a net zero magnetic helicity is supplied to the system by injecting equal and opposite values 
on sufficiently separated intermediate and small scales, the large scale magnetic helicity 
grows kinematically (independent of the magnetic Reynolds number) to approximately 
equal the sign and magnitude of the helicity injected on the intermediate scale. The 
oppositely signed magnetic helicity injected on the small scale subsequently drives a kinetic 
helicity driven dynamo through the secondary amplification of kinetic helicity. For the 
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driven case, in which the magnetic hehcity injection is continuous, hi grows beyond its 
kinematic value to a late time saturation value proportional to Rm for a closed or infinite 
system. For the free relaxation case, in which the driving is turned off, the secondary 
dynamo growth amplifies hi only slightly beyond its kinematic growth value of /i2(0)/2, 
to a maximum value that is independent of Rm and < h2{0). A real astrophysical corona 
likely never reaches the Rm dependent regime as Rm is extremely large in practice. The 
magnitude of the large scale field growth taken from the end of the kinematic regime is 
then the most relevant. For bihelical injection of equal and opposite signs on the k2 and ks 
scales, both injected with maximally helical fields, the ratio of magnetic energy on the large 
scale to that of the intermediate scale at the end of the kinematic regime is approximately 
ki/k2. 

When the appropriate initial conditions are used, the solutions of the generalized 
equations also reproduce the results of the previously studied simpler cases of unihelical 
magnetic relaxation^^ and the kinetic helicity driven dynamo. 

In the context of magnetically dominated astrophysical coronae of stars or disks, 
bihelical magnetic helicity injection is natural because such coronae are fed with fields 
by velocity driven dynamos operating in their rotators below. Nonlinear helical velocity 
driven dynamos produce helical fields of opposite signs on different scales^^' which can 
then emerge as bihelical fields into the corona. (Shear, not considered here, may effect the 
relative fraction injected on each scale^^' ^^.) The "large" scale and "small" scale of the 
rotator correspond to the "intermediate" and "small" scales of the corona respectively. 
The corona is therefore a natural environment for bihelical relaxation. The results herein 
suggest that helical large scale coronal fields can rapidly grow even when the coronae are 
injected with a net zero magnetic helicity. This is favorable for the in situ production of 
fields mediating coronal holes or MHD outfiows from stars and disks. 

Much more work is needed to incorporate the basic concepts and idealized solutions 
developed herein into realistic models of coronae. Shear, non-helical magnetic energy 
dynamics, and realistic boundary conditions should all be included. 
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Figure 1: Driven bihelical magnetic relaxation for ci — 0. Here h2 — —h^ — 1, ior all 

time, h = 1, k2 = 5, h = 20 and hi{t = 0) = 0.01. (a) hi, where Rj^ = Ry = 1000 
(dashed) Rm = Rv = 6000 (solid). The saturation value of hi depends linearly on 
Rm- (b) Same as (a) but for early times. The kinematic regime ends where the two 
curves split. Note that just after the kinematic regime the Rm = 6000 curve is below 
that of Rm = 1000 in contrast to the late-time regime in (a), (c) /i2 (short dashed, 
straight line) and h^ (long dashed). The curves for Rm — 6000 and Rm — 1000 
indistinguishably overlap, (d) h^ (straight dashed line) and h^ (oscillating solid line). 
Again neghgible variation for Rm — 1000 and Rm — 6000 cases, (e) Dimensionless 
EMFs: Q (solid humped curve), q (long dashed), and q2 (short dashed) for Rm — 1000. 
(f) hi when h^ << 1 is artificially enforced for Rm = 1000, and multiplying the q2 
term in (34) by 1 (dashed) and 0.5 (solid) to show the relative effect on hi. 

Figure 2: Same as Fig. (la) but with different values of the hi loss coefficient ci, 
and hi{0). (a) Ci = 0.1, hi{0) = 0.01 (b) Ci = 1, hi{0) = 0.01. These late time curves 
are insensitive to the choice of hi{0). 

Figure 3: Kinetic helicity driven dynamo. Here the system is driven such that 
hX = 1 for all times, and a seed /ii(0) = 0.01. These curves show that Eqs. (31-38) 
reproduce the results of Ref. 21 for the appropriate initial conditions, (a) Late time 
solution of hi for Rm = Rv = 1000 (dashed) and Rm = Rv = 6000 (solid), (b) early 
time solution for hi. 

Figure 4: Unihelical magnetic relaxation. Here the system is driven with h2 = I 
and /i3 = for all times, with a seed /ii(0) = 0.01. These curves show that the 
generalized set of equations reproduces unihelical relaxation of of Ref. 13. (a) hi for 
Rm — Rv — 1000 (dashed) and Rm — Rv — 6000 (sohd). (b) same as (a) but at early 
times (c) hi for Rm — Rv — 1000 (solid) and Rm — lORv — 1000 (dashed) (d) same 
as (c) for early times. 

Figure 5: Free bihehcal magnetic relaxation. Here hi{0) — 0.01, /i2(0) = — /i3(0) = 1 
but unhke the case of Fig. 1, there is no subsequent driving, (a) hi, where 
Rj^ ^ Ry = 2000 (dotted) Rm = Rv = 10000 (dashed) and Rm = 10000 and 
Rv = 2000 (solid). Note the slower decay for larger Rm and Ry. (b) Same as (a) 
for early times. Note the insensitivity to Ry at these times since the dotted and 
solid curves defined in (a) overlap, (c) /?.2 for Rm = Rv = 2000 (thin solid); h^ for 
Ryj = Ry = 2000 (short-dashed); /12 for Rm = Rv = 10000 (thick sohd); h^ for 
Rm = Rv — 10000 (long-dashed). Note that /i2 initially decays and /i^ rises and 
then the curves decay together. The solutions for /i2 and /i^ are insensitive to Ry 
for a given Rm. Note that /i2, h^ actually deplete faster for larger Rm (see (e) and 
text), (d) hs and h^ . At early times (barely visible), |/i3| decays and increases 
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and then the curves nearly overlap in decay: and for Rm — Ry — 2000 (top); 
/13 and hX for R^ = Ry = 10000 (bottom) /13 and for Rm = 10000, Ry = 2000 
(middle), (e) Blow-up of (0?) to show the offset between and with the curves 
are the solid curves on the right of each pair, and h-^ is the left dashed curve in each 
pair. The pairs from left to right are /is and h]^ for Rm = Rv = 2000; and for 
Rj^ = Ry = 10000; /i3 and for Rm = 10000, Ry = 2000. The offset accounts for 
the weak driven dynamo that grows hi slightly beyond its /ii = 1/2 kinematic 
growth value, (f) Dimensionless EMFs: q for Rm — Rv — 10000 (thick solid), q 
for Rm ^ Rv ^ 2000 (long dashed), ?2 for Rm ^ Rv ^ 10000 (thin sohd), ?2 for 
Rj^^Ry^ 2000 (short dashed). 
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